Quantum estimation of a damping constant 



Hannah VenzEl and Matthias Freyberger 
Institut fiir Quantenphysik, Universitdt Ulm, 89069 Ulm, Germany 

We discuss an interferometric approach to the estimation of quantum mechanical damping. We 
study specific classes of entangled and separable probe states consisting of superpositions of coherent 
states. Based on the assumption of limited quantum resources we show that entanglement improves 
the estimation of an unknown damping constant. 

PACS numbers: 03.67.Mn, 42.50.Dv, 03.65.Yz 



i> : 
o , 
o ■ 

Oh. 

<\ 

a^ : 



> 

a^ 

CD ■ 
cn 
o : 
i> ■ 
o ■ 

a ■ 

a • 

^ : 



I. INTRODUCTION 

Nonclassical quantum states can be used in various 
ways to tune up measurement resolutions beyond the 
standard limits [l|, [3] ■ The application of squeezed states 
of light in high-precision interferometry often serves as 
a prominent example [1, 0, Here the improved 

precision results from quantum noise reduction due to 
a specific superposition of photon number states. Hence 
appropriately prepared quantum superpositions offer one 
way to beat standard quantum limits. However, the 
quantum mechanical state space allows us to consider 
more sophisticated, nonclassical resources. In fact one 
can exploit the entanglement advantage of coupled quan- 
tum systems. The corresponding nonclassical correla- 
tions have been shown to enhance the measurement pre- 
cision Examples can be found in interferometry 
in quantum ellipsometry [§| , in atomic state spectroscopy 
, in the improvement of frequency standards [ll!, EM' 
and in optical lithography (isj . 

All these examples of quantum metrology are influ- 
enced by decoherence effects. In general any coupling to 
an uncontrollable environment will diminish nonclassi- 
cal enhancements. Therefore, it seems that decoherence 
has only a negative side when it comes to measurements 
based on quantum systems. However, counterexamples 
are known. Specific schemes for the generation of non- 
classical correlations require decoherence [ij]- Moreover, 
the fragility of quantum coherence effects can be used ef- 
ficiently to determine absorptions in optical sensor tech- 
nology [l5| . 

In the present contribution we shall describe an inter- 
ferometric estimation process for a damping constant. If 
we choose nonclassical input states for the interferometer, 
the estimation quality will be considerably higher than 
in the corresponding classical case. Moreover, we shall 
discuss how much of the effect solely relies on quantum 
superposition and what can be associated to quantum 
entanglement. 

The paper is organized as follows. We first introduce 
the model and the considered interferometric estimation 




FIG. 1: The setup consists of a two-mode input state trans- 
formed by a 50:50 beam splitter. We assume unknown 
damping only in the first mode with the damping constant 
7. A measurement apparatus (e.g., eight-port interferometer 
01 [13) registers eigenvalues of the EPR observables X and 
P, Eq. dH). 



scheme. Based on this scheme we shall then motivate 
different classes of nonclassical and classical probe states 
and discuss their properties in some detail. The different 
estimation results will be compared in order to determine 
the role of quantum superposition and entanglement for 
the quality of the described estimation scenarios. 



II. MODEL 

The schematic setup is shown in Fig. [T] It consists of 
two spatially separated modes A: = 1, 2 of the electromag- 
netic field, described by the mode annihilation (creation) 
operators ak {a\). They are connected by a lossless 50:50 
beam splitter. The corresponding unitary transformation 



C/bs = e" 



(1) 
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in general entangles a two-mode input state 

\^') = \^),®\^)^. (2) 

The degree of entanglement will depend on their spe- 
cial form. Behind this beam splitter, damping occurs 
in the first mode with the unknown damping constant 
7. Optical homodyne techniques can be used to measure 
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quadrature observables 



(3) 



for each mode k with Ok G [0, 2tt). In particular, eight- 
port homodyning [H,|T7[ can be appUed to determine the 
Einstein, Podolsky, and Rosen (EPR) observables [l^ 

X= -^(Xi-.T2), P=^{pi+P2), (4) 

where Xk = ifc(O) and pk = Xk{^). 

Our aim is now to show that such an interferometric 
setup allows us to characterize the damping. In par- 
ticular, we shall find that appropriate nonclassical input 
states, Eq. ([2]), which can be entangled by the beam split- 
ter, will enhance the quality of this characterization [l5|. 



A. Damping 

A common way to describe damping of an electromag- 
netic field mode is to consider its coupling to a reservoir of 
many modes, e.g., the radiation field of the vacuum. Here 
we will assume the usual linear coupling of the reservoir 
and the mode operators which models amplitude damp- 
ing. 

Behind the beam splitter this leads to the master equa- 
tion 



p(t) {~iu:\a\ak., +lCip, 



(5) 



k=l 



for a two-mode density operator p, under Born-Markov 
approximation in the zero temperature limit, where uj is 
the frequency of both modes. The dissipative evolution 
of mode 1 is described by the Lindblad term 



l^iP=l^ [^aipa[ - a\aip 



pa\(Xi 



(6) 



where 7 is the damping constant. Since we are consid- 
ering optical modes, the zero temperature assumption is 
reasonable. 



For an initial density matrix element 
in a coherent state basis one finds [22, 
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for the solution of Eq. ([5]) in the interaction picture. Here 
we have introduced the scaled damping constant k = jt. 



which allows us to compare the fidelity of the estimation 
using different input states, Eq. In particular we 

can choose input states which lead to nonclassical corre- 
lations between the modes of the interferometer. 

As mentioned before, we assume damping only in mode 
1, or in other words, the damping in mode 2 is supposed 
to be much weaker than the damping in mode 1. The 
second mode is thus basically considered as a perfect ref- 
erence mode. However, the more general case of indepen- 
dent damping in both modes can be analyzed analogously 
and leads to very similar results. To derive an equation 
for the average error obtained when estimating the un- 
known scaled damping constant k, we assume that we 
perform N measurements of an observable O, where Q is 
in general a function of the EPR observables X and P, 
Eq. ([¥]). The N measurements are supposed to be inde- 
pendent and carried out using identically prepared probe 
states, Eq. ([2]). The outcome of a single measurement run 
therefore consists of N eigenvalues {ili, = {^i} 
of the observable Ct. 

In order to find a simple expression for the estimated 
parameter k and under the assumption of weak coupling 
and short damping times, we linearize the expectation 
value 



d{n){K) 



(8) 



which contains k via the time evolved density operator 
p{k). We can then write linear function of {p}j{K) 



Co + ci(17)(k). 



(9) 



where the coefficients Cm follow directly from Eq. ([5]). 
Based on the N measurement results 57^ of a single run 
we will estimate the scaled damping constant 



Kest({^^j}) = Co + Ci 



1 ^ 
N ^ 



(10) 



As long as the linearization is valid this choice guarantees 
that /vest({^^i}) ~^ ^ fo'^ ^ ~^ We note, however, that 
this simplified approach can be refined. One possibility 
would be to obtain Hcstii^i}) from the optimization of 
an appropriate likelihood function [23] . In order to bring 
out the nonclassical features of our estimation scheme it 
is sufficient to work with the simple linearization. 



A. Average estimation error 



III. ESTIMATION OF AN UNKNOWN 
DAMPING CONSTANT 

After this discussion of our model, we shall now intro- 
duce the estimation scheme. We then derive an equation 



Considering identical experiments, the estimated val- 
ues Kost({f^i})i Eq. pO)) . will still vary from run to run 
according to a distribution which depends on the input 
states and has the expectation value k. If we want to 
compare certain estimation scenarios based on different 
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input states, we thus have to find a quantity defining the 
fidelity of a scenario. A simple possibility is the average 
width 



Ak^ = (k - Kcst({^^J)) , 



(11) 



as a function of the parameters defining our input states. 
The averaging is performed over infinitely many runs. 
Note that, for a single run, the quantity (k — Kost({f^i})) 
depends on the special set {fli, ^^n} of measured val- 
ues; see Eq. (fTO| . 

As mentioned above, the individual measurements in 
each run are independent and carried out using identi- 
cally prepared states. We can therefore write 

Ak2= j[] J dfl,Win,,K)] (K-Kest({f^J))' (12) 



does not involve any further resources to entangle input 
states. We therefore consider two scenarios comparable 
if the average total photon number in one run, 

ntot=N{n), (17) 

with the average photon number in one measurement, 

(n) = (*| (ni + na) |*) = {a\^i) + (18) 

is the same. We do not, however, fix how many mea- 
surements are performed in one run or, in other words, 
how many photons are contained in one measurement 
of f2. Using this idea of limited resources, the average 
quadratic error, Eq. (fT5|) . 



"tot \ an 



(19) 



with the probability distribution 

w{n,K) = {n\p{K,) \n) 



(13) 



given by the eigenstates of the observable Ct and the 
two- mode density operator p{k). Wc note 



dn W{n, k) = (fi")(K), (14) 



which is consistent with KestH^i}) = k. 

The average estimation error, Eq. (fTTjl . then reads 



A.2 = | [{Cl'){n)-{n)\n)) 



N ^ ' 



dn 



(15) 



where we have introduced the variance 



(16) 



and used Eqs. ^ and 

Considering, for example, an estimation based on the 
observable = X, suitable probe states should have min- 
imal noise in the two-mode X quadrature and the initial 
steepness [d{X) (k) / dn] \ ^=0 of the first moment should be 
large. 



B. Limited resources 

Since we want to compare the behavior of different es- 
timation scenarios we have to define which scenarios are 
considered comparable. In order to do so, we first assume 
that the special form of the input states Eq. 
has no influence on the preparation complexity. Further- 
more, the beam splitter is a passive element and hence 



has to be minimized as a function of the parameters de- 
scribing the probe states for a given number ritot- 



IV. PROBE STATES 

Our aim is to compare classical and nonclassical esti- 
mation scenarios based on the error given by Eq. (jl9p . 
We are therefore interested in finding appropriate ob- 
servables Cl, which are in general functions of X and P, 
Eq. and suitable input states, Eq. Guided by the 
expression for the average error we have to select those 
input states which show minimal noise in the correspond- 
ing observable. 

In order to work out the role of entanglement for the 
estimation of an unknown damping constant we regard 
both entangled and separable superposition states and 
compare them to what we will refer to as classical probe 
states. We will therefore discuss our setup, see Fig. [U 
with and without the beam splitter respectively and de- 
fine suitable probe states for each case. 



A. Entangled probe states 

To motivate the first set of probe states, we first con- 
sider our setup, see Fig. [U with the beam splitter in 
place but without any damping. So, in general, we will 
entangle the input states. For the identification of the 
parameters controlling the noise Afl^, Eq. (flB)) . we have 
to examine the probability distribution 



WiX,P)=\(X,P\UBs\'^) 



(20) 



which is determined by the input states transformed by 
the beam splitter and by the simultaneous eigenstates 
\X,P) of the commuting observables X and P, Eq. ([4]). 
We note that those eigenstates can be written as 



\X,P)=Ubs \p = P),(g>\x^-X), 



(21) 
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using the quadrature eigenstates \p)-^ and \x)2 of pi and 
X2- For the separable input state j^*), Eq. ((J), we then 
find 



W{X,P) = \{p = P\^)^\ |(x = ~X|V)2 



(22) 



for the probability distribution of the observables X and 
P, Eq. Note that each input mode only affects the 
distribution of one of the considered observables. The 
noise can be controlled independently. 

As shown before, for the estimation we have to find 
states that allow us to minimize the noise in the mea- 
sured observable 17. From Eq. we can see that 
ideal input states should consist of two quadrature states. 
Since quadrature states have infinite energy, i.e., average 
photon number, we rely on a more physical class of in- 
put states which can approximate those states to some 
extent. Clearly the first choice would be appropriately 
squeezed states. However, it has been shown, see, e.g. 
[24| . that by simply superposing coherent states, noise 
can be suppressed in single- mode quadratures Xk{Ok), 
Eq. ([3]). When we finally add damping, this has the ad- 
vantage that a simple analytic solution to the two-mode 
master equation, Eq. dS]), is known; see Eq. ([7]). 

Moreover, we restrict our analysis to the case of mea- 
suring the observable Cl — X, Eq. ([4]), only, since the 
principle can already be obtained from this reduced prob- 
lem. We thus conveniently consider superpositions of co- 
herent states in mode 2 and vacuum noise in mode 1 . The 
first class of input states will therefore be of the form 



a 
I — 
2 



a 

—I — 
2 



(23) 



properly normalized with the constant Ma and a G M. 
For those states, the probability distribution, Eq. ([22|) . 
has the simple form 

W{X, P) - e-^'-^' (l + cos W2aX]) . (24) 



It can then be shown in a straightforward way that the 
variance AX^ is given by [2^ 



2(1 



V2)- 



(25) 



See Fig. [2] The minimum, min {AX^} = 0.22, is below 
the classical value of 0.5 and is reached for \a\ = 1.6. The 
corresponding variance in P stays constant, AP^ = 0.5. 

For the states given by Eq. we further find that 
the distribution W{X,P), Eq. ([241), is localized at the 
origin of X-P phase space; that is, l^X^ = l^P) — 0. 
Since we want to see a strong influence of the damping 
on first moments, see Eq. ([T5|) . we have to displace this 
distribution in X-P phase space. We note the relation 



(X,P|i^l(Xo + lPo) 



(X-Xo,P-Po|, (26) 



where the phase T is a function of X, P, Xg, and Pq and 
£)fc(5) = e'^'^^-'^*'^'- (27) 
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FIG. 2: Variance AX^ Eq. as a function of the param- 

eter a. The fluctuations are determined by the input mode 2 
only. 



is the single-mode displacement operator of mode k. 
Moreover, the beam splitter transformation yields 



i?i(Xo + zPo) C/bs 



Xq + ^Po 
V2 



(28) 



Applying this combined displacement transformation to 
an arbitrary input state, Eq. ([21), will therefore not 
change the form of the corresponding probability dis- 
tribution W{X,P), Eq. but simply shift it to the 
two-mode phase-space point 



(X) = Xo, (P) = Po 



(29) 



Since we are only considering a measurement of the 
observable X, Eq. we displace the states according 
to 



(30) 



In other words, we increase the intensity in both input 
modes by the same amount. These input states are then 
transformed by the beam splitter. We will therefore use 
the states j^i) entangled by the beam splitter as suitable 
probe states for the estimation of the unknown damping. 



B. Separable probe states 

In order to distinguish the role of entanglement from 
the role of quantum superposition for the estimation of 
unknown damping, we now consider our setup, see Fig.[Tl 
without the beam splitter. Similar to the procedure in 
the preceding paragraph, we will define suitable separa- 
ble probe states, guided by the interferometric principle. 
Fig. [1] without damping. The probability distribution 



W{X,P) EE I (X,P I*) 



(31) 



is in this case directly given by the input state \^>), 
Eq. (21), and will not factorize like Eq. (j^ . In prin- 
ciple, the ideal minimal noise states are now given by 
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FIG. 3: Variances AP^ Eq. ((Ml) (dashed line) and AX^ 
Eq. (|33p (solid line) as a function of the parameter a. The 
dotted line represents vacuum noise. It can easily be seen 
that reduced noise in X strongly affects the width in P. 



eigenstates of the observables X and P, Eq. Q, as can 
be seen from Eq. ([51]) . Since those states cannot be pre- 
pared experimentally and above all are not separable, we 
again try to approximate them with superpositions of co- 
herent states. The restriction to a measurement of the 
observable fl = X, Eq. leads us to the symmetric 
input states 



l*ii)=Ar2 



a 
i — 
2 



2/2 



a \ 

-'2)2. 
(32) 



with the normalization Af^ and (^X") — (P) = 0. The 
form of the input states l^*!!) can be understood by 
noting that X, Eq. ([4]), consists of a sum of one- mode 
quadratures Xk- Therefore, noise has to be suppressed 
now in both input quadratures. 
The variances are given by 

1 ^,2 

AX^ = 



2 (l + e"V2) 



and 



AP^ = - 



V2 



2 (l + e"V2)' 



(33) 



(34) 



respectively. States of the form l^'ii), Eq. ((32)) . therefore 
allow for reduced noise in the observable X, Eq. ((3|), to 
the same extent as the entangled probe states Ubs I^i), 
Eq. ((2S|), discussed in the preceding section; see Eq. ((55|). 
However, as expected, reducing the noise in X leads to 
a broadening of the probability distribution, Eq. (|3ip. in 
P; see Fig. (3] The minimal value of AX^ = 0.22, for 
example, corresponds to AP^ = 1.5 which lies above the 
classical value of 0.5. The restriction to separable states 
leads to the feature that only one of the observables X 
or P, Eq. ((H, can be chosen to show reduced noise at a 
time [H. 

In analogy to the discussion in the preceding paragraph 
we can define 



as the final displaced class of separable probe states with 
(X) = Xo and (P) = 0. 

The two classes of probe states j^i), Eq. ((30p. with 
beam splitter and |$ii), Eq. ((55)) . without beam sphtter, 
respectively, will make it possible to analyze separately 
the influence of quantum superposition and entanglement 
on the quality of an estimation scenario in our scheme. 

Finally, we also compare these two classes of states 
to the most classical case with just one coherent state 
in each of the two input modes. Those states can be 
obtained by just setting a = in Eqs. ((50)) and ((55)) . 

Depending on whether the interferometer is used with 
or without beam splitter this leads to 



$111) = 



V2/1 



-Xr 



V2 



and 



|$iv> - |Xo)i® 10)2 



(36) 



(37) 



respectively. Before we come to the actual estimation, 
we shall now analyze the important features of our probe 
states in a bit more detail. In particular we shall outline 
how they are influenced by damping. 



C. Average photon numbers 

We first note the average photon numbers in one mea- 
surement, Eq. ((T8)) . for the different classes of input 
states. We obtain 



(ni) = (ni +7^2) |$i 



^0 + 4« 



^ tanh 



a 
T 



(38) 



for the input states |$i), Eq. ((30)) . and 



1 

(nil) = ^0 + tanh — 



(39) 



for the probe states |$ii), Eq. ((55)) . respectively. To gen- 
erate the same amount of noise reduction more photons 
are needed in the separable than in the entangled case. 
The number of photons is identical for both classes of 
classical states, (nm) = (niv) = X^. 



D. Influence of damping 

For estimating an unknown scaled damping constant, 
the effect of the damping on the probe states of course 
plays the crucial role for our considerations. According 
to Eq. p5p . two factors are of major importance: The 
variance AX'^{k), Eq. ((TO)) , and the first moment (X)(k), 
Eq. ([8]). For each of the considered probe states the first 
moment is given by 



|$ii) =i3i(Xo)|*ii) 



(35) 



(1)(k) =Xo 



-k/2 



(40) 
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FIG. 4: Variance AXj{k), Eq. (gT}, of the entangled probe 
states (solid line) and the variance AXfi{K), Eq. (|42p . of the 
separable superposition state (dashed line) as a function of the 
scaled damping constant k. Here the parameter a is chosen 
such that the states initially show maximal noise reduction. 



which basically shows the dissipative part of the damping 
mechanism. For the evolution of the variance, in which 
the decoherence is expressed, we find 



and 



AX2(«) 



1 _ a2(i + e-"/2)^ 

2 " 8(l + e°V2) 



(l+e--^) 



4(1 



V2)' 



(41) 



(42) 



respectively; see Fig. |4l It is clearly visible that the en- 
tangled state Ubs Eq. pO|). loses the noise reduction 
faster than the separable superposition |<i>ii), Eq. ([55)1 . 
However, due to the fact that we are considering asym- 
metric damping, i.e., only damping in the first mode, also 
the final values of the variances are not equal. Clearly, 
the variances of the classical coherent states stay con- 
stant, i.e., AXj2jj(k) = AXfy{K) = 0.5. 



V. ESTIMATION RESULTS 

We shall now compare the different estimation scenar- 
ios based on the basic assumptions described in Sec. IIIII 
We are thus searching numerically for the minimum of 
Ak^, Eq. p^ . under the constraint of a constant pho- 
ton number ntot 7 Eq. (fT7|) , for the various discussed probe 
states. In order to compare the estimation results of non- 
classical and classical input states, we define the relative 
improvement 



A4i 



(43) 



Here di {5u) denotes the relative improvement of states 
of the form I7bs Eq. dSO]), (|$ii), Eq. ^) compared 
to classical input states |$iii), Eq. The compari- 

son to classical input states I'&iv)! Eq. (|37p . is of course 
equivalent. 
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FIG. 5: Relative improvement Eq. (|43|l . in the esti- 

mation when comparing nonclassical probe states to classi- 
cal probe states. The dashed line represents Su, the com- 
parison of separable superposition states |"l?ii), Eq. psp . to 
classical input states |$iii), Eq. (|36p . and the solid line repre- 
sents Si, the comparison of entangled superpositions Ubs 
Eq. (|3Up . to classical input states. 



The obtained numerical results are shown in Fig. O 
for up to ritot = 20, using a scaled damping constant 
of K = 0.01. It is clearly visible that, independent of the 
number of available photons, applying either separable or 
entangled quantum superpositions leads to a higher qual- 
ity of the estimation. By using nonclassical probe states 
for the estimation, the error can be reduced by more 
than 50% when compared to classical input states. In 
addition, for our class of entangled superposition states 
Ubs Eq. ((30)) . the error is smaller than for the sep- 
arable superpositions Eq. ([55)1 . 

For just coherent states in each mode the error Ak^ 
scales like n^J^. However, for both separable and entan- 
gled superpositions the scaling depends on the number 
of available photons. Only from about ntot > 5 on the 
behavior is again governed by Ak^ cx n^Jt, which can 
be seen in the saturation of i5i and Su in Fig. [5l This 
effect can be explained by noting that for a fixed pho- 
ton number the average quadratic error Ak^ , Eq. ([19]) , is 
influenced by both the variance AX'^{k) and the initial 
steepness of the first moment [d{X){K)/dK]\^. For small 

photon numbers the suppression of noise in X, Eq. ([?]), 
has not yet reached its maximal value. For larger photon 
numbers all the additional photons are then, like for clas- 
sical states, being used to increase the intensity and thus 
the initial steepness of the first moment [d{X){K)/dK]\^. 
As discussed in Sec. llVCi less photons are needed to gen- 
erate the same suppression of noise when using entangled 
instead of separable superposition states. This leads to 
the smaller average quadratic error for entangled states. 

It is also interesting to note that as soon as the possi- 
bility of noise reduction is used it is on average best to 
use all available photons in one measurement. This can 
again be understood by recalling that a given fraction of 
photons is always needed for the suppression of noise. 

We conclude by stating that, in our interferometric 
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scheme, applying two-mode quantum superpositions im- 
proves the estimation of an unknown damping constant. 
Entanghng those superpositions by a beam spHtter even 
further improves this estimation. However, the main im- 
provement when compared to an estimation with classical 
input states is given by the reduced noise in a quantum 
superposition state. 

We further note that our estimation scheme can be 
considerably refined. Considering, for example, squeezed 
states or different observables could lead to an additional 
improvement of the estimation. By applying an adaptive 
estimation scheme [2^ one could possibly obtain a differ- 



ent scaling of the error Ak^, Eq. ([T9)) . with the average 
total photon number ntotj Eq. (jl7p . even for larger pho- 
ton numbers. 
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